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Abstract 

This work examines the effect of multiple photon emission on the quan- 
tum mechanical state of an electron emitting synchrotrotron radiation and 
on the intensity of that radiation. Calculations are done with a variant 
of perturbation theory based on the use of extended coherent states. A 
general formula is derived for the number of emitted photons, which al- 
lows for taking into account their mutual interaction. A model problem is 
used to demonstrate the absence of the infrared catastrophe in the modified 
perturbation theory. Finally, the elecron density matrix is calculated, and 
the analysis of this matrix makes it possible to conclude that the degree of 
the electron's spatial localization increases with the passage of time if the 
electron is being accelerated. 
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1. INTRODUCTION 

The effect of radiation on the path of charged particles in a synchrotron has 
already been analysed (see, e.g., Ref.l and the literature cited therein). The 
analysis is based on the classical Lorenz-Dirac equation or on the solution of 
a kinetic equation whose coefficients are the probabilities of quantum transi- 
tions between various stationary states of an electron moving in the magnetic 
field of the synchrotron. Calculations have shown that in the absence of fo- 
cusing in the magnetic field there is increase in the radial fluctuations of the 
electron path and an increase in longitudinal fluctuations of the electron's 
momentum with the passage of time. In a focusing magnetic field, in the ini- 
tial stages of electron acceleration, the presence of radiation leads to radiative 
damping, which damps radial and vertical oscillations the so-called radiative 
damping effect). Lately, research are focused on the analysis of equation of 
Lorenz- Dirac type in problem not necessarily related to synchrotron radia- 
tion (see, e.g., Refs. 2 — 4). New solutions of Lorenz- Dirac equations have 
been found for some special cases, and additional arguments from quantum 
electrodynamics are invoked to eliminate nonphysical solutions. 

Despite the indisputable value of the results obtained by solving Lorenz- 
Dirac equations, it must be noted that some important properties of the 
states of a particle emitting radiation, properties that are not directly related 
to the path and do not directly influence the above effects — like an increase 
in radial fluctuations of the path - - are excluded from these results. This 
is true, in particular, of the evolution of the particle's wave packet, which 
affects the radiation and hence the radiative friction and the path. 

The present work is an attempt to use a modified perturbation- theory 
approach to examine the effect of multiple photon emission on the evolution 
of the wave packet of a particle, in particular when the particle emits syn- 
chrotron radiation. Most papers devoted to the quantum mechanical theory 
of the synchrotron radiation ignore this aspect. It is usually assumed that 
the particle emitting radiation has a wave function given by the solution of 
the Dirac equation. However, if the emitted radiation is taken into account, 
the particle is only a part of of the quantum mechanical system and its state 
cannot be described with the completeness that is possible in principle in 
quantum theory. 

A common approach to describing the states of particles that are members 
of a large system is to use the concept of the density matrix. This work 
demonstrates that the evolution of the density matrix suggests that a particle 



2 



goes with the passage of time into states that are more and more localized, 
with the particle motion described by the laws of classical mechanics with 
ever-increasing accuracy. Thus we have additional support for the validity 
of using the Lorenz-Dirac equation along with a clearer undestanding of the 
incomleteness of the physical picture described by this equation. 

The problem of the structure of the wave packets of emitting particles is 
related to the classical model of a distributed electron studied by Lorenz (see, 
e.g., Ref. 5). In quantum electrodynamics this model leads to the well-known 
problem of ultraviolet divergence, encountered in the calculations of the mass, 
charge, and energy of an elementary particle. Furthermore, the renormaliza- 
tion of charge in quantum electrodynamics reveals the internal inconsistency 
of the traditional Feynman formulation of perturbation theory (see, e.g., Ref. 
6). It would be useful to follow the changes in the difficulties encountered 
by classical electrodynamics initiated by changes in the perturbation theory, 
to establish which of the above problems is invariant, so to say. It might 
turn out that in the modified theory some of these problems can be resolved 
without resorting to additional hypotheses. This might then lead to a new 
direction in the development of quantum electrodynamics and the theory of 
quantized fields in general. The presen work uses a model to show that at 
least in relation to the infrared catastrophe, the adopted modification in the 
theory does not lead to problems characteristic of the traditional form of the 
theory. Calculations are based on a general formula that describes the mu- 
tual interaction of the emitted photons as a manifestation of the nonlinearity 
inherent in quantum electrodynamics. 

2. EMISSION OF PHOTONS BY A CLASSICAL CHARGED 
PARTICLE 

We write the Hamiltonian describing the interaction of a free electromag- 
netic field and a particle carrying an electric charge Z (here we use atomic 
units: Ti = 1 and |e| = 1): 



where the current density vector j is a function of coordinates and time. The 
vector potentual operator is specified in a three-dimensional transverse gauge, 




c 



(i) 



a,q 



qr 



(2) 
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as a standard linear form in the creation and annihilation operators (fl q 
and f aci for photons in states with polarization a(a = 1,2), momentum q, 
and energy u = cq. The polarization vectors e aq have unit length and are 
orthogonal to q. The coupling constants g q = (2ttc 2 /cjQ) 1 / 2 contain the 
normalization volume Q, which does not enter into the final expressions and 
thus can be put equal to unity. 

The total Hamiltonian is the sum of the free photon Hamiltonian 

#0 = E W /Iq/aq + COUSt. 
aq 

and the Hamiltonian (1). We pass to the interaction picture for field opera- 
tors, for example, 

A(t) = e lflat Ae- l%t . 

The equation describing the evolution of the wave vector of the photon field, 
\i) in the interaction picture, 

if t \t) = B int (t)\t), (3) 

has in the given case an exact solution in the form of the direct product of 
photon coherent states, 

\t) = II expHXa,q - /aqQaq + /ctqQaq] l*o), (4) 
a,q 

where the initial state vector coinsides, to within the arbitrary phase factor, 
with the vacuum state of the photon field: 



| to) = e l ^°\vac), 0o = const 



to 

t 

XaS) = J MQ«q(0<&q(0]< 

to 

and j q (t) is the Fourier transform of the current density. 

Using the exact solution (4), we can calculate all quantities of interest. 
For instance, the mean number of photons created by time t is given by 

= l««o(*)l 2 - ( 5 ) 
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Note that this formula yields the mean number of emitted photons only as 
t — > oo, since creation of a photon requires a time interval c/q long, which 
tends to infinity as q — > 0. In what folows we use the interpretation of n aq . 

Suppose that the photons are emitted by a point particle carrying electric 
charge Z and moving along a path r = ro(t). Then 

jqW = Zvo(t)exp[-iqr Q (t)}, 

where vq = ro(t). After summing over over polarizations we can reduce the 
time derivative of the number of photons (as t —> — oo) to the form 

T . £ nffiW = -ff, 2 / drlMt ~ \r\/2 + r/2)v„(t - \r\/2 - r/2)- 

a—i,z _ oc 

i(qv (* - \r\/2 + r/2))(qv (* - |r|/2 - r/2))]* 

exp[iuT - iq{r {t - \t\/2 + r/2) - r (* - |r|/2 - r/2))]. (6) 
Applying this equation to the case of synchrotron radiation, we obtain 

d_ »< 2 00 



v 2 00 1 
^n^ = Z 2 -|^ / dr[coso;or - -cos 2 6(cosu;or + cos((2t - |r|)o;o))]* 



-00 



exp{icjr — 2iqR cos© sin^- cos(uo(t — |r|/2))}, (7) 

where uq = eH /7rac, t>o = i?wo,H is the magnetic field strength, R is the 
orbit's radius, m is the particle mass, and 7 is the Lorenz factor. The angle 
6 is the inclination of the vector q to the orbital plane. 

The expression (7) is periodic in time, with period To = 2ti/uq. Averaging 
over one period, we obtain 

j 2 00 rl / \ 

I E n { Sl = Z 2V \g\ j dr e™ ± cos 2 6 J 2 2gi? cos 6 sin ^f] + 
dt a c J 1 2 \ 2/ 



—00 



cos wor ( 1 — - cos 2 G^ Jq f 2gi? cos G sin 



(8) 



Next, we allow for the fact that for any periodic function F(r), 

OO T 

j e^ T F(r)dr = £ e ma;T ° / e wr F(r)dr, (9) 



-00 n=-oo 
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where the sum of exponentials can be transformed into a sum of delta func- 
tions: 

OO Orrr OO 

E e inwT ° = — E S(u-n'u ). 

n=-oo -L n'=—oo 

Combining this with (9), we can transform (8) to the following form: 

— J2n [ ai = Z 2 ^g 2 q £ - / dxe 2mx {cos2x J {2qR cos 6 sinrc)(l--cos 2 6)+ 

fit a C n=0 7T J Q 2 

^ cos 2 6 J2(2qR cos 6 smx)}27i5(u — tiuq). (10) 

Since u = cq > 0, the sum in (10) goes from to oo. 
Next we have 7 

J e 2lfiX J 2 v(2a sinx)dx = ne^^J^^a) J v+fi (a), 
o 

and the recurrence formulas 

2n 

Jn+l(z) + Jn-l(z) = J n (z); 

Jn-i(z) = -J n (z) + J' n (z); 

71 

Jn+l{z) = ~Jn{z) ~ J' n (z). 

25 

As a result, Eq. (10) becomes 



^£ng> = Z 2 g 2 q t 27r6(u - nu )[tg 2 O J n (-cos6) + cos 6)]. 

(it a n=0 C C C 

(11) 

Equation (11) can be used, in particular, to obtain the well- known Schott 
formula. Thus, for the mean intensity of synchrotron radiation the semiclassi- 
cal theory yields results that coinside with classical results. The semiclassical 
theory provides additional information (in comparison to that provided by 
classical electrodynamics) only in the sence that it makes it possible to cal- 
culate the fluctiations in the number of the emitted photons, their mean 
energy, and total momentum. Similar results can be obtained for the case in 
which a charged particle moves along an arbitrary path 8 . In all cases we at 
classical formulas for the mean intensity of radiation emitted by the particle. 
Moreover, calculations of the mean electromagnetic field that accompanishes 
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a charged particle moving in vacuo lead to well-known expressionsfor the re- 
tarded potentials 9 . This agreement between the semiclassical and classical 
theories forms the basis for a more accurate quantum mechanical theory of 
interaction of radiation and an emitting particle. 

3. QUANTUM MECHANICAL THEORY 

We consider the interaction of an electron and the radiation emitted by 
that electron. We pass to the furry representation and write the wave oper- 
ator of the electron in the form of an expansion in the stationary states of 
type (A2) (see the Appendix): 

where we have excluded the antiparticle operators, since allowing for the 
contribution of particle-antiparticle intermediate states leads only to small 
corrections to the phenomena considered. The creation and annihilation op- 
erators, (i[ and must obey the standard Fermi commutation relations. 
The current density operator is approximately (without allowing for electron- 
positron pair contributions) given by 

j a (t) = cfta a $ = cE^^^l^Hr)^^"^'^ (12) 

ft' 

(from now on a, b\ = 1, 2, 3 label the projections of vectors on the Cartesian 
coordinate axes). 

We construct the operator 

^\t)=Zr ( q ,t)e- i ^)p(0)^ (13) 

where i*o(q, t) is a vector (which needs to be determined) that depends on the 
momentum transfer q and time t, and is the "zeroth" density operator 
at time t = 0: 

Pq 0) = E dlJk+w 

k,cr 

Here o\ a and d^ are the creation and annihilation operators for an electron 
in a state with momentum k and a projection of the electron spin on the z 
axis that takes the values a = ±|. Note that the operator (13) is selected in 
a form that satisfies the charge conservation law. 

We require that the running mean Fourier transform of the operator (12) 
coincide with the expectation value of (13): 

(^ q WW=r (q,t)e-^o(q.*)( t |^o)| t) _ (M) 
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Then i"o(q, t) must be approximately equal to the mean position of the particle 
at time t. We define the deviation of the current from the "zeroth" value to 
be 

Aj q (t)=j q (t)-j(°)(t). 

This deviation will be used to build the interaction operator in the new 
representation. The above transformation is convenient because the operators 
(13) commute at different times: 

[j q 0,a W,jT(*')]=0. (15) 

Using (13), we write the electromagnetic interaction operator as a sum of 
two terms, H int (t) = H^(t) + H^(t), where 

H^(t) = --Jf°\t)A(t)dV; (16) 

= -- J m)A(t)dV. (17) 

c j 

Then, by virtue of (15), the equation 

= ff»»(t)|t) (18) 

has an exact solution in the form of a direct product of the vectors of extended 
(or modified, in terminology of Ref. 10) coherent states, 

\t) = II eX P ("% q - UQU + /aq^aq) |0), (19) 
a,q 

where (at to = 0) 

QaS)=i-}dtK^\t l )e iuJt \ (20) 
c o 

XaS) = --J {SlqCOQaqM - 5UO<LqWK- (21) 


The initial state vector |0) is the direct product of the vacuum state of the 
electromagnetic field, \vac), and the vector of the initial state of the moving 
particle, \(f>), described by the wave function (f)(r),i.e., |0) = \(f),vac). 

We have chosen to = to be zero rather than — oo due to the fact that, as 
further calculation show, the temporal sequence of changes in the state of a 
moving particle that interacts with the field of the radiation it emits is highly 
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important. In this approach there are sure to be problems assosiatad with the 
inteaction turning on, which violates the charge conservation, and with the 
generation of virtual radiation, which is the consequence of such violation. 
To avoid the need to discard fictitious terms, one can resort to turning the 
interaction on slowly by replacing the constant Z with a slowly increasing 
charge Z(l — e~ et ), where e is small. The charge buildup time Ti n = e _1 must 
be long compared to uj~ , but short compared to the observation time (here t 
must be much longer than u ). After we establish a method for evaluating 
the integrals for some definite value of q, we can extend it to any other value 
of q. 

If we ignore the corrections generated by H^t, Eq. (19) fully solves 
the problem of calculating the physical quantities of interest. In particu- 
lar, instead of (5) we have 

»aq(*) = (*l^q(*)/aq(*)l*) > 

which at ro(q, t) = ro(£) leads to a result coinciding with (5). Thus, if we 
ignore the current variant of the quantum mechanical theory differs 

from the semiclassical one in calculations of the mean number of the emitted 
photons or the energy of these photons only when ro(q, t) differs from ro(t). 

At the same time, corrections due to H^ t , can be obtained for any con- 
venient choise of the vectors ro(q, t). Let us put ro(q, t) = ro(t). Then 
Qaq(t) = Qaq(t)p q , where Q a q(t) is specified by its semiclassical expression 

Z 1 

QaS) = J dt'e* a(l vo(t') exp{iujt' - *qr (O}- (22) 
c o 

To construct a new "modified" perturbation theory in H^ t , we introduce the 
zero-order evolution operator 

U (t) = exp{J2 Qaq(t)fl q ~ Q ] aq{t)faq ~ %qW}- 
a,q 

Then (19) can be written as \t) = Uo(t)\Q). We also introduce a new repre- 
sentation of operators: 

A(t) = Ui(t)A(t)U (t). (23) 
The state vector \t > in this representation obeys the equation 

d\t >= H^m > ■ (24) 
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Allowing for (24), we can reduce the expression for the mean number of 
photons to 



x 1 i I tn—l 

n aq (t) = n^(t) + £ (-i) n J dti J dt 2 ... J dt n * 

71=1 



<0|[.4/ii(*)/a^ (25) 

4. CALCULATING CORRECTIONS IN THE MODIFIED THE- 
ORY 

Writting the series in Eq. (25) explicitly, we find that the expansion 
contains terms proportional to even powers of Z. We collect the leading 
terms of this type, which contain Z 2 as a pre-exponential factor. Such terms 
exist in n^(t) and in the first and second of the series (25). Note that here 
Z is also contained in the exponents entering into Uq and Uq. We calculate 
the first commutator in (25) via the following auxiliary formulas: 

faS) = (/aq + QaS)) ^ /JqW = (/lq + QlM ^ 

[f aq ,uo(t)} = uo(t)Q aq (ty, \JUMt)] = -ukmU^ 

IfL Uo(t)} = U (t)QUt); [/ aq , Ui(t)] = -W (t)Q aq (t). (26) 



We put 
so that 



Baq(t) = e aq Aj_ q (f), 
H$(t) = -- E 9 9 (f aq B aq (t)e-^ + P m Bl q (t)e^) . (27) 



a,q 



Using (26), we can perform th following transformation: 

[/i q (*)/«qW > ^W] = 
t%(tl)l{Fa« + @c«(t,t$ (28) 

where Q a q(Mi) = Q aq {t) — Q aq it\). Since the operators Q aq {t) already 
contain Z as a factor, the leading terms emerge as a result of the commutation 
of the photon operators and H^ t : 

IfU^faS)^^)} - |^t7 t (t 1 )(e-^5 aq (t 1 )(/ aq + Q aq (t,t 1 ))- 
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^(/Iq + QUMllKttl))^!). (29) 

If we average (29) over the initial state of the system by employing the equal- 
ities 

/aq&0(tl)|0)=Qaq(tl)&0(*l)|0); 

(0\ul(t 1 )fl q = (0\ukti)QU t i)^ 

we get 

(0\[P aq (t)f a S),W t (ti)}\0) = ^ 9 MUkti){BaSi)QaS)e- iLJtl - 

^a q W5i q (ti)}[/o(ti)|0). (30) 

In calculating the next corrections in (25) we emmediately discard terms 
that contain pre-exponential factors with Z raized to a power greater than 
two. This means that when we plug such terms into the second and subse- 
quent terms of the sum in (25) into the expression for the first-order commu- 
tator, we can immediately discard terms containing the operators Q aq and 
Qt. In the resulting expressions, the operators / aq and /t can be freely 
interchanged with the operators Uq and Uq, since their commutators contain 
heigher-order corrections in Z, which we have just discarded. 

In view of this, all terms in which the annihilation operators f aq are to 
the right of other /- operators, or in which the creation operators /t are to 
the left of other /-operators, must be dropped. In the remaining terms the 
operator products faqfa'q' m ust be replaced by the commutators 5 aa >A(q — 
q'). By performing these transformations we reduce the leading terms that 
appear when we write the double commutator on the right hand side of 
Eq.(25) explicitly to the form 

-f^ g 2 (cft(*i)e-^ 

^o t (^)5 aq (^)t7o(t2)^o t (^i)5i q (^i)^o(^i)e^ (<1 ^ 2) )^ (31) 
Collecting all terms of the specified order,we get 

Z 1 

naS) = \Q*S)\ 2 - i-9qj dh^Ul^ie-^B^Qa^t)- 

c o 

c o o 
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5 aq (t 1 )[7o(^i)^o t (^)5i q (^)[7o(^)+ 

Ukt2)Bi q (t 2 )Uo(t 2 M(^ (32) 

If we now write all terms in (32) that appear because of plugging the explicit 
expressions for B aq (t) into (32), collect like terms, and do the necessary 
canceling, we arrive at the final result: 

naS) = —9 2 q j dti J dt 2 (0\ul(ti)* 

c 

eaJU^i)^o(ti)t/o t (^)e;J q (t 2 )t/o(t2)|0)e-^-^). (33) 

Note that at deriving (33), we did not take advantage of the fact that ro(q, t) 
is independent of q, with the result that the formula still holds in the general 
case, in which ro(t) is replaced by ro(q, t) in(22). 

5. NUMBER OF PHOTONS 

We assume that in the expansion of the initial state vector of the particle, 
|0), the expancion coefficients Cj^ in states with definite momentum |kj) have 
a peak at ko, and decrease as kj ko devates from ko, by the Gauss law 

C ki = (2^i) 1 / 2 (2 7 r^) 1 / 4 exp[-^ ( 5f/4 - (k l± - k 0± ) 2 ^/4], 

where k = (k±,pi), ko_L is time- dependent, and 5i and 5± - are longitudinal 
and transversal packet widths (relative to the z axis). This representation 
follows from the study of electron states in a magnetic field in the Appendix. 
For relativistic electrons, the momentum uncertainty in the initial state is 
much less than the momentum proper. In real calculations of the numbers of 
emitted photons via (33), it is preferable to represent the current operators 
as expansions in states with definite momentum at a given moment in time, 
with a time dependence characteristic of plane waves. In the present paper, 
this approximation is justified by the fact that due to the strong effect of 
the radiation on the particle's state in the comoving reference frame, an 
effect exeeding the one produced by the external field, we can ignore the 
quantization of levels in the time dependence of the operators. Indeed, even 
the classical theory of synchrotron radiation predicts that the mean energy 
of the photons emitted by a particle is much greater than luq. In view of 
this, the mean difference in particle energies before and after photon emission 
proves to be much greater than the separation between the levels of transverse 
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motion. Under these conditions, allowance for level quantization in the time 
dependence of the operator can only lead to small corrections of order 1/n 
(where n ~ 7 3 - is the mean ratio of the frequency of the emitted photon to 

As a result oof the action of electron operators, the vectors k and ki in 
the current operators in (33) are transformed into the vectors kj — Aq , where 
Aq = Eq s , with qi, q2, ... the momenta of emitted photons. Replacing the 

s 

given expression with kj(t) = k— Ak(t), where Ak(t) is the mean momentum 
loss by the particle by the time t, and plugging it into all the cofactors in 
(33) that are not in the exponential, we get 

n aq (t) = %9\\ dh J dhe-^-V £ |C k J 2 (e aq v^(q, *i))« q v^(q, f 2 ))* 
c o o k i^' 

(vac; ki,ff|C^(ti)4(*l)%(tl)C^(t2)^q(*2)%(*2)|k 

u a;vac), (34) 

where 

Pq(t) = E 4-q,a'i,aexp{i(£ k _ q - £ k )t}, 
k,a,cr' 



o a o h k\{t) + 



£i + mc 2 
\ £ ■ + mc 2 



o\k\(t) - q b ))a a 



w 



ui 



with £i = £\ il (t), z'i = £ ki (t)-q, and the summation over repeated indices is 
implied. The term coresponding to a' ^ a describes emission processes 
accompanied by electron spin flip. Further simplification is possible if the 
exponents in the density operators in (34) are transformed according to 



(35) 



where /i(q, q'J are unspecified functions. In this approach, different photons 
are assumed to be almost independent, since otherwise we would have to 
speak of a strong correlations between the emission of two separate photons, 
which agrees neither with the semiclassical theory nor with the calculations 
below. In an approximation that is linear in Aq, for q <C ki we have 

Mi(q, q'J ~ (V£ ki -q - V£ ki )q' s « -qq^/m^, (36) 
where 7^ = £^ t /mc 2 . As q' s — > 00, the function //(qq'J ceases to depend on q' s . 
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Using the methods of calculating means employed in Ref.9, we get 

t t 

n aq {t) = JdhJ dt 2 E |aj 2 QL q (^i,Ofeq(^,Oexp[-P, q (^i^2)], (37) 



where 



Z 1 

Qiaqit^') = i~9q [ e* q v i)(7 / (q, t')exp[iut' - zqr i>CT ,(q, t')]dt', 



c 

u 

with rj ;C r'(q, t) = Vj ;Cr >(q, t). The exponent in (37) is given by 

Pid(t h t 2 ) = E \Qif3q'{h^')\ 2 (1 - exp[-^(q,q%]) + 

\Qiw{t 2 ,(j')\ 2 (1 - exp[i/^(q,q%]) - 
QW^ 1 ' ^OQwfe </) (1 - exp[-i/^(q, q%]) * 

(l-exp[i^(q,q%])]. (38) 

Obviously, 

\imP iq (t h t 2 )^0; 
lim P iq (*i,* 2 ) -> 0. 

Equation (37) contains the desired corrections to the semiclassical expres- 
sion for the number of emitted photons. It assumes its semiclassical form 
for |-Pi q (ti, £2)1 <C 1. From a physical standpoint, this difference betwwen 
the formulas is due to the fact that in (37) we allow for interaction of the 
emitted photons, while in the semiclassical theory this factor is ignored. The 
probability distribution for the number of emitted photons in each state does 
not obey the Poisson law any longer, which a reflection of the nonlinearity of 
electromagnetic phenomena in the quantum theory. 

Obviously, an equation like (37) can be used to study arbitrary motion of 
a particle, not just an electron in a synchrotron. To do so, we merely redefine 
the quantities Vj CT /(q, t), which in the simplex case can be approximately 
calculated for the mean of the vector and averajed over spin (in this case, 
the velocities v(q, t) and the function (38) no longer depend on the indices % 
and a'). 
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Let us estimate P q for the case in which the velocity v(q, t) is constant 
and equal to vq: 

Q aq (t) = -g q ^^e^-^)\ 
c u — qv 

Plugging this into (38), we obtain an expression that is logaritmically diver- 
gent, due to the slow decrease in the integrands as q' — > oo. This fact is the 
manifestation of ultraviolet divergence, often encountered in electrodynam- 
ics. In contrast to Feinman's perturbation theory, ultraviolet divergence does 
not lead to a catastrophe: it only means that (in contrast to the predictions 
of the semiclassical theory) a uniformly moving particleis not accompanied by 
transverse photons. This example is a clear demonstration of the dependence 
of the way in which the ultraviolet divergence depends on the perturbation 
theory employed. A detailed study of this problem lies outside the scope of 
the present paper, where we use the standard method of introducing a cut- 
off momentum q c ~ mc to remove the ultraviolet singularity. The resulting 
expression for the absolute value of the function (38) proves to be small and 
varies very slowly (logaritmically) with t\ and t^- An explicit estimate of the 
function (38) for v(q, t) constant will be made in the next section. 

6. INFRARED ASYMPTOTIC BEHAVIOUR OF THE NUM- 
BER OF PHOTONS 

Let us consider the asymptotic behaviour of the function (37) as u = qc — > 
0. In classical electrodynamics (see, e.g., Ref. 11) and in the semiclassical 
theory there is a characteristic frequency dependence of n aq as u — > 0, namely, 
n aq ~ 1/aA Hence, upon integration with respect to momenta, the total 
number of emitted photons diverges logaritmically at the lower limit. Will 
allowing for the effect of emission on a state of the emitting particle (as in 
Eq. (37)) influence this pattern? To ansver this question, we examing a 
model problem in which a charged particle moves at constant velocity Vi 
and, colliding at time t% > with a point scatterer, suddenly changes its own 
velocity by a small quantity Av = V2 — Vi, |Av| <C Vi, and then proceeds 
to move at constant velocity V2. The requirement that this jump in velocity 
be small simplifies all calculations considerably. Moreover, since a jump 
in velocity implies infinite acceleration, various nonphysical effects are to 
be expected. The requirement that the velocity jump be small makes the 
velocity almost a continuous function, so that such effects can be ignored. 
When Eq. (37) is employed in calculations, there is the problem of the 



15 



interaction suddenly turning on at the initial moment in time, which violates 
charge conservation, and of generation of fictitious radiation, which is the 
consequence of such violation. To avoid the need to discard fictitious terms, 
one can use the procedure developed in Sec. 3 to turn the interaction on 
slowly. 

Since the particle is assumed to have a definite velocity, we drop the 
subscript i in (38) and replace the vectors v ia with the current value of 
the velocity. We calculate the resulting functions P q (ti,t 2 ), assuming that 
Tin C ^3 < To this end, we first estimate the quantities f(ti : t 2 ) = 
expfizqq'ti^/rn/y] on the right-hand side of Eq. (38). Since a photon is 
emitted when the particle changes velocity, we consider the neighbourhood 
of the point t\ = £3, t 2 = £3, assuming that £3 ~ jom/q 2 . The vector q' is 
the momentum transferred from the moving particle to the emitted quanta 
of electromagnetic field (photons). The mean value of this momentum is of 
order ra|Av|, so that at q <C ra|Av| the ratio q'/q can be large. Thus, the 
absolute value of the exponent in f{t\,t 2 ) in the range of parameters under 
investigation is large, and the exponentials are rapidly varying functions that 
make a negligible contribution to (38). Eliminating these contributions from 
the outset, we reduce (38) to the simpler form 

P^hM) = T,[\Q^(ti)\ 2 + \Q^(h)\ 2 - 

QwWQMfo) (1 + e-^ 1 -^)]. (39) 

We now calculate the function (39) explicitly for v(q, t) = Vo = const. In this 
case, assuming that q c is much less than the mean momentum of the emitting 
particle, we calculate the integral with respect to q' and obtain 

UhM) = ^/^ ^^ (iSiMt! - **))+ 

% Si((<j 2 + - t 2 )) + 2C- Ci(<j 2 |ti 

Ci(|cj2 + - t 2 \) + In (u 2 \u 2 + u^h - £ 2 ) 2 )), (40) 

where n' = q'/V, u\ = g c n'q/m7, U2 = (c — nVo)g c , Si(£) and Ci(^) - are 
the sine and cosine integrals, and C = 0.5772 • • • - is Euler's constant. The 
function (40) vanishes at t\ = t 2 and slowly increases with the time difference 
At = \ti — t 2 \. In the nonrelativistic limit at large At ^> l/cq c , the function 
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(40) can be arrpoximated by the expression 



Pq(*l,*2) 



2Z 2 v 2 

37TC 3 



in sign(^ - t 2 ) + 2 (C + ln(cg c ) + ln|*i - * 2 Q] • (41) 



We remark on the smallness of the coefficient of the expression in square 
brackets. As At increases, the real part of (41) increases logarithmically, but 
the characteristic buildup time proves to be exponentially large, so that the 
function (41) can be considered small over the entire range of its arguments. 

Now let us estimate the number of photons emitted by the electron in the 
entire course of its motion for the nonrelativistic case. Integrating by parts, 
we find, for instance, that 



naqif) = ~ l ^9 2 q / dtie atl v(ti)e 



— w^ti+iqr^i). 



u - qv(t 2 ) + idP q (t 1 ,t 2 )/dt 2 



exp{icut 2 - iqr(t 2 ) - P q (ti,t 2 )} 



t 2 =t 
t 2 =o~ 



J dt 2 exp{iut 2 - iqr(t 2 ) - P q (ti,t 2 )} 



d 



J ac\ 



V(*2) 



dt 2 \u - qy(t 2 ) + idP q {t l M)/dt 2i 



Allowance for the value of the first term inside the square brackets at the 
lower limit is unjustified because of the violation of charge conservation at 
t — > 0. If we turn the interaction on slowly, then this contribution is zero. 
First we integrate by parts with respect to ti, using the same ideas that we 
used in integrating with respect to t 2 . We obtain 

Z 2 



n 



aq 



e aq v 2 



u - qv 2 + idP q (t,t 2 )/dt 2 



t 2 =t 



+ 



g\ J dt 2 J dti exp{iu(t 2 - h) + iq(r(ti) - r(* 2 )) - P q (t h t 2 )}* 



d 



eaqV(t L 



d 



dt\ \-u — qv(ti) — idPq/dti dt 2 \u — qv(t 2 ) + idP q /dt 2 



(42) 



where we have discarded the rapidly oscillating terms, which contribute noth- 
ing to the overall expression for the number of emitted photons. The first 
term on the right-hand side of Eq.(42) corresponds to the part of the trans- 
verse field that follows the moving particle, and is related neither to change 
in the particle's velocity nor to the radiation. Hence in all calculations of 
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the characteristics of the radiation that follow, we allow only for the second 
(integral) term. 

In calculating the time derivatives in (42) we encounter continuous and 
delta-function terms, with the latter being a reflection of the discontinuity in 
velocity, the derivatives Q aq and P q . For instance, 



d 



V(*2) 



dt 2 \u - qv(t 2 ) + idP q (t 1 ,t 2 )/dt 2 



i®(h-t 2 ) 



d 2 P q (h,t 2 ) 



e; q vi 



dt\ (uj-qv 1 +idP q (tht2)/dt 2 )' 



iS(t 2 -t 3 ) 

s(t 2 - 1 3 ) 



d 2 P q (h,t 2 ) 



< q V 2 



dtl 



{u-qv 2 + idP q {t 1 ,t 2 )/dt 2 ) 2 

-1- 

>v 2 



u - qv 2 + idP q (ti,t 2 )/dt 2 



t 2 =t 3 +o 



aq 



(43) 



t2=t 3 -0 



u - qvi + idP q (ti,t 2 )/dt 2 

Here 0(0 is the Heaviside step function. 

The relationship between the continuous and delta-function terms in (43) 
can be evaluated as follows. The total contribution of the ©-functions can 
again be calculated by parts, which again results in a delta-function contri- 
bution multiplied by the magnitude of the discontinuity of the integrand at 
t 2 = ts. This jump includes the second derivative of P q as a factor whose 
order of magnitude can be estimated to be the product of the first derivative 
and the mean value of the frequency of the emitted photon. The latter cannot 
exceed the energy lost by the moving particle, and it is therefore proportional 
to the small parameter A = viAv/^. Clearly, allowing for the continuous 
terms in (43) would mean allowing for the next terms in the series expansion 
of the integrals in A. the leading term is still the contribution of the delta 
function, the only contribution we consider. 

Using the condition that the interactionis turned on slowly, we find that 



QaS) = —9q 



e aq V l 



u — qv x 

For ts < t the result is different: 



-e^-^i)*, Tin « t < t 3 . 



Qaq(t) = —g q 



e« q vi 



e aq V 2 



u — qv x u — qv 2 



e i(w-q Vl )i 3 _^ 



<q V 2 
( 

lo - qv 2 



J(u>-qv 2 )t 
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For Ti n <C t 2 < ^3 we have 



dU 



= -*£<%^i) -^(e; q ,v 1 )e i ^-^^(l+ 

/3q' LC 



qq 

exp{-iqq (ti - t 2 )/m^}) H Qp q >{t 2 ) exp{-iqq (ti - t 2 )/m'y} 

777.7 



Finally, for £2 > ^3 we have 
^2 W L •" 



*2-*3- 



= E -2^-<#sin{(u/ - q V 2 )^-^}(e^v 2 )(-: 



e/Jq'Vi 



e/3g^v 2 
u/ - q'v 2 



qq 

7 — 

7777 



a;' — qVi 



*Ba'(ti)Q0d'(t 2 ) exp{-iqq'(*i - t 2 )/m^}- 



^J q /(^)-^(e^v 2 )e^-'^' V2 ) <2 (l + exp{-7qq / (7 t 1 - * 2 )/m 7 }) 

Note that dP^ti,^) / dt 2 is continuous in 74. 

Let us calculate the delta-function contribution to the integrals with re- 
spect to t 2 in (42), letting t — > 00 : 



(00) = — ^ I efei exp[-7o;(7 t i - * 3 ) + 7q(r (7 t i) - r (* 3 )) - P q (*i, £3)]* 



2 00 
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aq 



5 



e aq v(ti) 



d*i La; - qv(*i) - idPq(ti,t 3 )/dti - qv 2 + idP q (ti,t 2 )/dt 2 



e aq V 2 



e aq V l 



u - qvi + idP q (ti,t 2 )/dt 2 

Now we integrate with respect to ti, again limiting ourselves to delta- function 
contributions. Allowing for the fact that P q (t 2 ,ti) = P*(ti,t 2 ), we obtain 



Z 5 



n aq {oc) = —g q 



e aq V 2 



u - qv 2 + idPqfo, t 2 )/dt 2 



t 2 =t 3 +o 



<qVl 2 



(44) 

w-qvi + iaP q (t 3 ,*2)/^2 

This equation solves the problem. If we neglect the derivatives of P q in the 
denominators, (44) coincides with the standard expression for the number 
of low-frequency photons emitted in a collision, the expression that can be 
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derived in classical electrodynamics 11 and in quantum electrodynamics if we 
use standard perturbation theory 6 . Let us estimate the derivatives of P q in 
the denominators in (44). We have 

Oh t2=t3+Q c q , g (a; -qvij 

where Oi(g) is of the first order in q . For q small, noting that 



t 2 =t 3 -o c q' Q -qvij 



at 2 

and that V2 ~ vi, we obtain 

Z 2 

n aq (oo) = 

where 



e aq V2 e aq vi 



u — qv2 + A a; — qvi + A 



(45) 



A = 2^^, ,y, XV f v (46) 

In the nonrelativistic limit v q <C c , from (46) we obtain A « 4Z 2 t> 2 (/ c /37rc 2 . 
Equation (45) does not contain the infrared singularity. A reviation from the 
n aq ~ l/o; 3 law with decreasing cj begins at an energy of order A. The lower 
the energy of relative motion of the charged particle and point scatterer, the 
lowe the aforementioned energy. This estimate also holds if the velocity of the 
particle changed not suddenly but over a time interval that is short compared 
to the time of production of a low energy photon. 

7. CALCULATING THE DENSITY MATRIX FOR THE CASE 
OF SYNCHROTRON RADIATION 

We now use the above approach to calculate the density martrix of an 
emitting particle. The exact expression for the density matrix in the repre- 
sentation realized by the transformation (23) has the form 



7 (x,x',t) =< tW(x,t)iJ>(x!,t)\t > . (47) 

We calculate (47) in the first approximation, replacing the vector \t > by 
the initial state vector |0 >. Using the Baker-Hausdorff rule with proper 
transformation of the evolution operators Uo(t) and Z7q (t) , and the operators 
Qaq(t) in the form (20), we easily find that 

7 (x,x , ,t) = 7 (x,x , ,t) e xp[-S(x-x',t)], (48) 
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where 70 (x, x',t) = 0*(x, t)0(x', t) is the value of the density matrix tthat 
does not account for emission and is determined by the wave function 0(x, t) 
of the exactly described state of the electron in an external magnetic field. 
The function S(x. — x',£) in the exponent is given by 

5(x - x',t) = £ |Q aq (*)| 2 [l - e^"')] (49) 

q,a 

which vanishes at x = x'. As |x — x'| — > 00, the function (49) acquires the 
maximum value, equal to the total number of photons emitted by the given 
moment in time. 

The mean momentum of the particle can be evaluated as follows: 

p(*) =Po + */ \^t)\ 2 VS(^x.',t)\ x=x/ d 3 x. (50) 

This means that the gradient V / S'(x, x', t) determines the rate of decrease of 
the mean particle momentum due to emission of photons. If the initial state 
was stationary, |0(x,£)| 2 does not depend on time. In this case, the mean 
force acting on the particle is 

F 6 = if |V>o(x,t)|VS(x,x',t)| x , =x d s x. (51) 

The calculation of the function 5(x, x', t) for the case of synchrotron radiation 
is similar to the calculation of the proton production rate in Sec. 4. Noting 
that S actually depends on the difference r = x — x', we obtain the value of 
s averaged over one period: 

72 



Z r 00 r 711)n 

S(r,t) = t— / d6sm6 £ nu cot 2 6 J 2 (— sin<9)+ 

r ' n=l ' '' 





■%J' 2 (-^^sin(9) [ 1 — Jo (r^^ sin 6$ sin 9) exp(ir^^ cos6*o cos<9) ) , (52) 
& c J V c c j 

where #0 is the polar angle of the vector r with respect to the axis perpen- 
dicular to the orbital plane, r = |r|. In the ultrarelativistic case the following 
approximate formula is more convenient: 



2 2 / 3 Z' 



2, , 00 



S(r,t) =t *~ I Zj U ° f d^ l/2> f d6 sin6* 



To 

1 

,2 



{ctg 2 # Ai%/2) 2 / 3 (l - ^ sin 2 0)1+ 

c 2 
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j ^2/3 Ai /2 [(,/2) 2 / 3 (l - -| sm 2 6)]}* 
^1 — Jo(sin6*o sin 6 —^) exp(icos#o cos 6 — , (53) 

where = e~ 3 , 0- = tt/2 - e, 6*+ = 7r/2 + e and 1/7 <C e <C 1; Ai(z) is the 
Airy function defined in the Ref.12 and Ai'(z) is its derivative. Obviously, 
the imaginary part of the averaged expression for S given by (52) and (53) 
is zero. 

According to Ref.10, the density matrix (48) describes an ensemble of pure 
states (in the sense of von Newmann), whose properties are determined by 
the behaviour of e~ s ( T,t \ The expansion of the matrix (48) in the density 
matrices of these pure states can be written 

7 (x,x',t) = / d 3 a$* a (x,t)$ a (x',t)N a (t), (54) 

where $ a are the wave functions of the pure states, and N a (t)d 3 a are the 
probabilities that these states are realized at the given momentum in time. 
The functions $a(x, 1) are proportional to the products of the wave function 
0x, t) and the wave functions x( x — a, ^) , where x( x , t) is the solution of the 
integral equation 

G(r, *) = e -' s( "'^ = J x*(x - a, t) X (x' - a, *) d 3 a. (55) 

But what about the existence and uniqueness of the solution of this equation? 
If we write (55) in the Fourier representation 

G-q(*)=x£(*kq(*), 

the absolute value of the desired function is uniquely defined, but not the 
phase. However, this uncertainty is a direct consequence of the translation 
invariance of Eq.(55), whose general solution, therefore, has the form 



Xq = e'VG-q (56) 

with arbitrary real a. A solution exists if the Fourier transform G q is a 
real nonnegative quantity. That it is real follows directly from the fact that 
lmS(r,t) is odd and Re5(r,t) even under insversion; the nonnegativity fol- 
lows from the fact that ReS'(r, t) increases monotonically with distance r. 

The effective size of the localization region for the initial state in the orbital 
plane is Sp ~ ^JR/^vq (see Appendix). The quantity Sp is usually much larger 
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than atomic dimensions. The localization region for the initial state along 
the magnetic field is infinitely large, which is due to the initial uncertainty in 
the z-component of the momentum. The latter is obviously determined by 
the macroscopic parameters of the actual experimental layout. 

Radiation can substantially alter the picture, and lead to spatial localiza- 
tion of the emitting particle in a region whose size is of the order of atomic 
dimensions. To estimate the rate of variation of the widths of the states <& a (t) 
with the passage of time, the function S(r, t) was calculated numerically for a 
set of parameters characteristic of the FIAN-60 synchrotron (E = 0.68GeV 
and R = 2m). 

8. MAIN CONCLUSIONS 

The perturbation theory developed in this paper has made it possible to 
establish that certain fundamental problems of quantum electrogynamics are 
not invariant when the type of perturbation theory is altered. With respect 
to infrared divergence, this is shown by emloying the simple example of an 
emitting particle that undergoes a sudden change in velocity. The results 
have been obtained for the nonrelativistic case, since the study of emitting 
relativistic particles requires a detailed analysis of the ultraviolet asymptotic 
behavior of the integrand in (38). 

The method of removing ultrfviolet divergences by introducing a cutoff 
momentum, which was adopted in the present paper, is not covariant under 
Lorenz transformation, and therefore cannot be used in a consistent rela- 
tivistic theory. But even preliminary studies show that in the new approach 
the problem of ultraviolet divergence is not catastrophic, in contrast to the 
case in ordinary perturbation theory of quantum electrodynamics. It is to 
be hoped that further research will lead to progress in understanding this 
problem. 

Density matrix calculation have shown that reduction of the spatial dimen- 
sions of the localization region for the emitting electrons to atomic dimensions 
can be achieved over a macroscopically long time interval r c of some tenths 
of a second. Can the present theory be applied to such long times? The sit- 
uation is complicated by the fact that in the course of one orbital revolution, 
the particle is subject to a solenoidal electric field that balanced the loss of 
energy to photon emission. If we assume that this field acts during a time 
interval so short that it only accelerates the particle's wave packet as a whole 
and is unable to change the particle's internal parameters substantially, then 
there is no reason why to do estimates we cannot extend the theory to the 
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entire duration of the particle's motion in the synchrotron. 

The time r c is much shorter than it takes the packet to spread due to the 
nonequidistant nature of the spectrum of the transverse-motion levels. What 
is observed is an anisotropy in the packet's width: the packet is most strongly 
squeezed perpendicular to the magnetic field, and least strongly parallel to 
the field. The considerable elongation of the packet in the direction of the 
magnetic field is obvious. 

The posibility of strong spatial localization of the emitting particles means 
that if the acceleration cycle in the synchrotron is long enough, the motion 
of the particle can be described to high accuracy by the equations of classical 
mechanics. Nevertheless, this does not mean that the intensity of the radi- 
ation must agree with the prediction of classical electrodynamics. Indeed, a 
localized state in quantum mechanics is completely different in its properties 
from a localized state in Newton's classical theory. The justification for using 
Newton's equations of motion to calculate the paths followed by wave packets 
is provided by Ehrenfest's theorem, but the decisive factor in calculating the 
intensity of the radiation is the momentum of the particle, rather than the 
position. In quantum mechanics, a state with a definite momentum is com- 
pletely delocalized, and in this way differs substantially from states of type 
$ a • There is thus no way in which we can intuitively interpret calculations of 
the characteristics of radiation using classical ideas. The characteristic com- 
mon feature of the formulas derived in this paper is the fact that allowing for 
the mutual interaction of the emitted photons reduces the radiative intensity. 
A similar result was obtained by Landau and Pomeranchuk 13 , who studied 
the radiation emitted by charged particles moving in continuous media (the 
Landau — Pomeranchuk effect). The physics of this phenomenon amount to 
the fact that random collisions of an emitting particle with particles of the 
medium can reduce the path length over which the radiative indensity builds 
up coherently. Something similar is observed when photons are emitted into 
vacuum: multiple emission of photons can mimic the multiple collisions in a 
continuous medium that lead to a reduction in radiative intensity. 

The numerical calculations were done together with V.A.Aleksandrov. 
The work was supported by a grant from the Russian Fund for Basic Re- 
search (Grant No. 97-02-16058). 

APPENDIX: COHERENT STATES OF A RELATIVISTIC 
ELECTRON IN A UNIFORM MAGNETIC FIELD 
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Let 

A = (-iyH , jrH ,0), 

where Hq is the strength of the magnetic field directed along the z axis. The 
motion of an electron in such a field, which obeys the Dirac equation, has 
been the topic of numerous studies (see, e.g., Refs. 14-16). The solution 
given below differs from the well-known one only in some details. 

We introduce the lowering operators for the two independent oscillators: 

51 = ^y/rnu L (x + ip x /mu L - iy + p y /mu L ); 

52 = ^mu L (x + ip x /mu L + iy - p v /mu L ), 

where lul = |e|Ho/2mc is the Larmor frequency. The frequency of the first 
independent oscillator is twice the Lormor frequency, u)\ = 2ul, while the 
frequency of the second oscillator is zero. The set of lowering and raising 
operators (which are conjugates of lowering operators) satisfies the standard 
Bose commutation relations. The operators a\ and a\ describe the orbital 
motion of an electron in a magnetic field, while the operators and a 2 
describe the position, fluctuations and other characteristics of the center of 
the osculating circular orbit, whose mean radius is R. 
We next introduce the matrix operator 



D = 



p z —2iy/mujL Si 

2iy/mu L a\ -p z 

The energies of the electron's quantum states are 



^ m 2 c 4 + p 2 c 2 + 4:ULmc 2 (ni + a + -), (^-1) 



where the label £ = (ni,ri2,cr,p) simply indicates the set of quantum num- 
bers in parentheses. The ni,ri2\ = 0, 1,2, ... label the quantum states of the 
independent oscillators, with n\ being the principal quantum number. We 
denote the projection of momentum on the z axis by p. The discrete variable 
a takes two values, ±1/2 , corresponding to two possible projections of spin 
on the direction of the magnetic field. The bispinor describing a stationary 
state of an electron in a magnetic field is given by 



1 ( JEe + mc 2 (fc(r) 

^ t] --wX%^^ ]e e ' (A2) 
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where v^( r ) is a spinor of the form 

^(r) = e^^==(a\r(4rifoAp)- 
Vni!n 2 ! 

Here 

(p Qa (p) = ^mu L /7T exp(-mu L p 2 /2)xa, 

with 

p 2 = x 2 + y\ Xi/2 = (1,0), x-i/2 = (0,l), 

and the normalization length along the z axis is taken equal to unity. 

An arbitrary solution of the Dirac equation is a linear combination of 
bispinors of type (A2). Just what linear combination corresponds to the 
initial state of an electron in the synchrotron? For standard values of syn- 
chrotron parameters ( for example, for the FIAN-60 synchrotron), the mean 
value of n\ is very large (of order 10 13 ), and if the assumption that photons 
are emitted largely independently of one another is true, so is Poisson's law. 
In this case, the expected value of the relative fluctuation of the number n\ 
is extremely small, A = Ani/ni ~ 1CT 6 . Hence, essentially all terms of the 
desired linear combination can be expanded in powers of A, with the result 
that 

E$ « E a p + (2u L /j (Tp )Ani, 

where E ap is the value of E% at £ = (ni, H2, &,p), and 7^ = E ap /mc 2 . We 
see that the spectrum is essentially uniformly spaced, with the leveles being 
separated by the mean orbital frequency u ap = 2^/7^ of the electron about 
the magnetic fiald. When the relative fluctuation of n\ is small, we can put 
E^ ~ E ap in all nonexponential factors. 

The linear combination corresponding to the above requirements has the 
form 

1 / \jE ap + mc 2 (p ap (r, t) \ 
if>*p(r,t) = -j=— p / f) exp(-iAE ap t), (AA) 

where AE Gp = E Gp — uj ap rTi, and 
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with a\(t) = a\exp(iu ap t), and a\{t) = a|, where ai and «2 are constant 
phases; the momentum along the z axis is assumed equal to p. The compo- 
nents of the current density vector in the state (A4) are 

2 C 2 

f ap = -^-y/mu; L ni\ip ap {r, t)\ 2 sm(u ap t + an) 
2c 2 

j v = --^VmJJtfh\ip ap (Y, t)\ 2 cos(u ap t + cti) (A.5) 



J ap 



is = --^w,t)l 2 . 

The packet's rms width in the radial direction in the state (A4) is deter- 
mined by the radial behavior of the function (fo a {p) and can be estimated to 



be Ap = y2Rc/ E ap . In the azimutal direction, the stationary states of type 
(A2) are completely delocalized. Indeed, in these states the angular momen- 
tum is well-defined, and by virtue of the uncertainty relation for action-angle 
variables, they cannot be localized in angle. 

In contrast, the state (A4) has no definite angular momentum, but its 
angular width is limited, and is of order A(f> ~ 1/ y/W\ in the azimutal direction 
(we assume that the uncertainty in the position of the orbit's center is much 
smaller than the orbit's radius, so that ri2 <C ni), which after being multiplied 
by the orbit's radius yields a distance roughly equal to Ap (for the FIAN-60 
synchrotron this distance is about one micrometer). 

The packet width along the z axis is governed by such macroscopic pa- 
rameters of the device as the diaphragm width, and for this reason it can 
exceed the radial or azimutal width many times over. In this case the packet 
can be represented by a linear combination of states of type (A4): 

i(r,t) = E^( r -i)- (4-6) 

p 

where the constants C p satisfy the normalization condition and guarantee, 
e.g., a Gausian dependence on the z projection of the momentum with mid- 
point at p = 0: 

C p = {2n5l) l ' A e-^\ 

If we assume that the spatial width of the packet along the z axes is much 
greater thanthe radial width, then ^> Ap, and in this case the state (A6) 
is associated with a small symmetric ellipsoid elongated in the direction of 
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the magnetic field and revolving in this orientation ina circular orbit about 
an axis parallel to z. To estimate the time of packet spread in the radial or 
azimutal direction, we must keep the next term in the expansion of the energy 
in powers of Ani. This yields the value of the time of packet spreading due 
to the nonequidistant levels of transverse motion, t\ ~ joR 2 / An\. Here 70 is 
the Lorenz factor for the electron beam in a synchrotron. For the FIAN-60 
synchrotron the time T\ was estimated to be about ten seconds. 

References 

[1] I.M. Ternov and V.V. Mikhailin, Synchrotron radiation [in Russian], En- 
ergoatomizdat, Moscow (1986). 

[2] A.I. Nikishov, Zh. Exp. Teor. Fiz. 110, 510 (1996) [JETP 83, 274 (1996)]. 

[3] V.I. Ritus, Zh. Eksp. Teor. Fiz. 110, 526 (1996) [JETP 83, 282 (1996)]. 

[4] G.F. Efremov, Zh. Exp. Teor. Fiz. 110, 1629 (1996) [JETP 83, 896 
(1996)]. 

[5] D.D Ivanenko and A. A. Sokolov, Klassische Feldtheorie, Akademie- Ver- 
lag, Berlin (1953). 

[6] V.B. Berestezkii, E.M. Lifshitz, and L.P. Pitaevskii, Quantum Electro- 
dynamics, 2nd ed.,Pergamon Press, Oxford (1982). 

[7] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Sums, Series and 
Products, Academic Press, New York (1980). 

[8] G.M. Filippov, Poverkhnost' No.3, 52 (1997). 

[9] G.M. Filippov, Poverkhnost' No. 12, 73 (1995). 

[10] G.M. Filippov, Zh. Exp. Teor. Fiz. 101, 1634 (1992) [Sov. Phys. JETP 
74, 871 (1992)]. 

[11] L.D. Landau and E.M. Lifshitz. The Classical Theory of Fields, 4th ed., 
Pergamon Press, Oxford (1975). 

[12] M. Abramovitz and I. Steguin, Handbook of Mathematical Functions, 
New York: Dover, 1965. 



28 



[13] L.D. Landau and I.Ya. Pomeranchuk, Dokl. Akad. Nauk SSSR 92 , 735 
(1953). 

[14] V.V. Dodonov, E.V. Kurmyshev, and V.I. Man'ko, in Proc. of the P.N. 
Lebedev Physics Institute [in Russian], Vol. 176, Moscow (1986), p. 128. 

[15] V.G.Bagrov, D.M. Gitman, and V.D.Skarzhinskii, in Proc. of the P.N. 
Lebedev Institute [in Russian], Vol. 176, Moscow (1986), p. 151. 

[16] V.G.Bagrov, D.M. Gitman, and V.D.Skarzhinskii, Preprint P.N. Lebedev 
Physics Institute No. 101, Moscow (1986). 



29 



